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\ Abstract. We quantify the conditions at which "composites", the resonances 

If} • and bound states d, He 3 are produced at RHIC. Using Hubble- like model for 

late stages, one can analytically solve the rate equations and also calculate the 
relevant optical depth factors. We calculate also the modification of p masse 
and width, and predict a radiacal shape change of a. 
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1. Introduction 



This talks is based on two papers, one with G.Brown [ ^ and one ongoing work 
with P.Kolb [|2|, who will report a part of it related with late time evolution here. 
fH \ Their common goal is to understand what are the conditions which determine the 

timing of production of observable resonances, and estimate at what conditions this 
happens. We have emphasized observable above because there are many resonances 
produced in the system but unobserved in the final invariant mass spectra because 
of re-scattering of their decay products. The development of simple analytic model 
of the kinetics of resonance production/absorption, as well as evaluation of the 
"optical depth" integral. 

An old but still important idea is "matter modification" of hadrons, providing 
the experimental test of the conditions in question. As such we discuss modifica- 
tion of the mass and shape of two classic ittt resonances, p and a. Recent STAR 
observations of these effects are reproduced (for p) and predicted (for a). 

The issue of production of nuclear fragments is also an old one, at RHIC reduced 
to d, He 3 and their antiparticles. In literature those are studied either by statistical 
or coalescence models, which left open many important issues. First of all, like 
resonances the observable fragments must escape all interactions with any particles 
in order to survive at the end. The second point is that this process is production- 
rate-limited, thus it can lead to non-thermal quasi-equilibria. Furthermore, new 
element is the consistent evaluation of the 2-to-l production rate itself in a recent 
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work by Ioffe et al [0]. 

2. The optical depth factor 

Let me start with the simplest pedagogical points about the observability of reso- 
nances and fragments, or "composites" as we may call them collectively, for brevity. 

In the next section we will discuss rate equations which can be solved and 
determine the number of composites N(t) at time t. The "observability condition" 
of a resonance can be written as 

v visi ue{t) = TN{t)exp (- j v(t')dA (1) 

where the the l.h.s. is the production rate of visible resonances, T is the resonance 
decay width and the exponent is the optical depth factor containing integrated v(t), 
the combined scattering rate for all decay products. The N(t) decreases with time 
due to expansion and cooling, while the exponent changes from at early time to 
1 at late times: so the product naturally has a maximum at the time t m such that 

1 ^ =0 (2) 



N(t m ) dt 

This condition means that for observable resonances the freezeout condition is dif- 
ferent from that for stable particles and reads: the rate of their number change 
is equal to the absorption rate of all the decay products. For example, for p and 
a we should not know their scattering rates but just that of two pions. Since for 
short-lived p and a the first factor is close to overall expansion rate of matter at late 
time which follows from Hubble-like late-time regime d\ogN(t m )/dt w 3/t, and the 
second is the same, we conclude that 'visible" p and a are produced at the same 
time. 

The formation rate for a fragment made of A nucleons is made by some coa- 
lescence, and such rate is obviously proportional to nucleon density to that power, 
~ (nN(t)) A . After it is produced, however, it still has very small probability to 
survive. Assuming that the destruction rate for a fragment va ~ Ai/n, where vn is 
a scattering rate for one nucleon, one finds that for A-fragment the time distribution 
is approximately the A-th power of the same universal function 



fragments (^) 



njv(t)exp(—J Vf[(t')d£^ 



(3) 



So, the maximum of production of any visible fragment happens at the same time for 
all A. Furthermore, the width of the distribution over production time decreases 
as A grows, as 1/ \J~A. 



Freezeout of Resonances 



3 



3. Solving the rate equations 



The equations themselves describing dynamics of resonances are well known, gener- 
ically they contain the sink (the decay) and the source terms 



dn(t, r) 

at 



-Tn(t, r) + S(rii 



(4) 



(where for expanding source the time should be understood as proper time in the 
rest frame of all volume elements.) In many papers in literature (e.g. [|3]) the source 
term is ignored citing "instantaneous hadronization" , but (especially for resonances 
we consider) it is not true: in fact the primary generation of resonances die out long 
before the "observable" ones are born. 

We use an approximate power fit of the source time dependence J d 3 rS = 
TN f^-J Its power can be related to fireball expansion. If the volume V(t) ~ 
l/n(t) ~ t a the integrated source is proportional to V(t)[nTr(t)] N where N is the 
multiplicity in resonance decay ( N = 2 for a, p) . The pion number is "chemically 
frozen", N n = V{t)n n (t) = constif) from which we conclude that the source term 
power is P — (N — l)a (for a, p and other binary resonances P = a). An example 
of such equation solved is shown in FigGJ 





Fig. 1. (a) The time dependence of the r/io-meson density, starting from chemical 
equilibrium at t = 5 and 10 fm. Plotting its ratio to the pion density in (b) one 
observes the transition from hydro to free streaming regimes and that the p density 
decreases power-like rather than exponentially, because of the source term. 



We then evaluate the optical factors for pions and nucleons, using realistic re- 
scattering rates, with chemically frozen composition, using papers by Hung and 
myself [ 5 and by Tomasik and Wiedemann Example for the final time dis- 
tributions for visible p and d is shown in Fig [3 Note that both distribution have 
maxima we discussed above, and that the "visible" d are indeed produced very 
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late. This is our main point: we are speaking about a very dilute matter, after the 
freezeout of all the basic ingredients of the fireball. 
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Fig. 2. The time dependence of the visible r/io-meson (a) and deuteron (b) 
production, for the r — point in central AuAu collisions at RHIC. 



4. The resonance modification 

It has been argued over the years that in matter the resonances should be modified, 
with shifted mass, increased width and even significantly changed shape. With 
the very late stages of RHIC collisions relevant, we can now access very dilute 
matter in which those effects must be calculable in the lowest order of the density, 
providing a benchmark test to all such discussions. In such case hadron modification 
is expressed in terms of their forward scattering amplitude My(i = 0, s). Note that 
the scattering amplitude is complex, and that this approach gives both the real 
and imaginary part of the dispersion law modification, also known as the optical 
potential. 

There are two major theoretical approaches to the issue discussed in literature, 
to be called an s — channel and a t— channel one. The former approach assumes that 
the scattering amplitude is dominated by s-channel resonances which are known to 
decay into the i + j channel. For most mesons such as n, u>, p, K in a gas made of 
pions such calculation has been made e.g. in related with such up resonances 
as ai or N p resonance TV* (1520) [0 for p at SPS. Note that the signs of the effects 
are opposite in those examples, as seen from the following table: 

The majority of the particles in the matter are Goldstone bosons n, K, r\ which 
do not interact at small momenta. However attraction between other particles is 
there. Using a simple expression for the mass shift one gets Sm 1 ^ « —28MeV due 
to all B + B. An additional shift 8m v p — WMeV comes from scalar exchanges 
between p and all other vector mesons p,uj,K*. The main difference between the 
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Table 1. A set of resonances considered 



Name/Mass 


Width 


Branching 


Mass Shift 


ai(1260) 


400 


0.6 


-19 


a 2 (1320) 


104 


0.7 


-15 


ifi(1270) 


90 


0.4 


+ 1.6 


K 2 (U30) 


100 


0.087 


-0.4 


N*(1520) 


370 


0.2 


10 



two mechanisms of the mass shift discussed above is that the t-channel attraction 
is not associated with the broadening, while the s-channel resonances increase the 
width by about 50 MeV. On the other hand, there is a "kinematic" effect working 
to the opposite direction. The negative mass shift discussed above automatically 
reduces the width, both because of the reduced phase space and also due to the 
power of p in the P-wave matrix clement. The magnitude of this effect for the 
predicted mass shift is 

A 777 

,5r p = 3-^r p » -50 MeV (5) 

m p 

So, inside the accuracy these two effects cancel each other. 

The invariant mass distributions in pp and mid-central AuAu of the 7r + 7r~ 
system, with a transverse momentum cut 0.2 < p t < 0.9 GeV have been measured 
by STAR see also the C.Markert's talk here. I would not have time here to 
discuss the shift in pp (see [H]), and I only comment that in AuAu the p peak is 
found to be shifted by additional ~ —40 MeV in mass, but the width is the same. 
This agrees well with estimates above. 

The same approach should of course be applied to many other resonances. 
For more narrow resonances, like K*, we expect smaller shifts, while for wider 
resonances like a we predicted a complete change of shape. At small freezeout T as 
100 MeV sigma was predicted to be deformed into a much more narrow structure 
at mass of about 400 MeV, see figures in [[I]. Exactly such a peak has been seen 
by STAR 1 . Another confirmation of very late freezeout and low T, the a/ p ratio 
strongly grows toward central collisions. We are waiting for quantitative analysis 
of these data with great interest. 

5. Fragment coalescence 

The issue of coalescence, such as p + n — > d, was discussed in many papers over 
the years, and authors struggled with the question how to calculate its rate. In 
particular, it is clear that when the level crosses zero the wave function at the 

1 l have seen it first at this workshop, on the day after my talk, shown by Gary Westfall in a 
plot of the balance function as a function of qinv with a peak he said "nobody understands". It 
was rather good timing between the prediction and its experimental confirmation. 
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origin vanishes, and so the production should do so too. And if the production rate 
is small compared to two other relevant rates, v a bs and dlogN/dt, there is never 
thermal equilibrium and one should not use statistical models. 

Significant progress has been made in recent paper by Ioffe et al [0| who have 
pointed out how to use consistently the in-matter widths of all particles and obtain 
the production rate. We are now incorporating it into the picture of expansion and 
the optical depths discussed above, and hope to get quantitative results for d, d 
spectra soon. 
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